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In metallic nanoparticles, the different electronic environ- 
ment seen by each magnetic nucleus produces a distribution 
of Knight shifts of the NMR frequencies which is observed 
as an inhomogeneously broadened lineshape. We study the 
fluctuations in the local density of states for s electrons at 
the Fermi energy in a simple LCAO model. We resort to a 
Matrix Continued Fractions calculation of the Green's func- 
tions. Results show that line broadens asymmetrically and its 
shift decreases as the particle size or the temperature dimin- 
ish satisfying a universal scaling function. However, for very 
small particles, surface states become relevant to determine 
a lineshape that departs from the universal scaling behavior. 
These trends are consistent with the observed tendencies in 
Cu and Pt particles. 
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I. INTRODUCTION 

Small clusters of atoms of|, metallic elements have 
interesting physical propertieso which originate in the 
strong interference effects that appear when the coher- 
ence length of the electronic excitations is larger thfui 
the size L of the particle. As already noted by KuboH, 
this would allow the manifestation of the geometry de- 
pendent quantization of electron energies. This condition 
leads to the very active field of mesoscopic physicsQ which 
involves a deep revision of traditional ways of calculation 
and interpretation of transport and spectroscopidj prop- 
erties. On the other hand the diverse technical applica- 
tions of metal particles as well as composites of metal 
particles has pushed the dejsjelopment of a wealth of sam- 
ple preparation techniquesECI which allow varied system 
characteristic. It has been recognizecfflD that NMR can 
be a useful tool to characterize metallic nanoparticles and 
that further insight in the comprehension of their prop- 
erties requires the development of models which allow a 
quantification of spectra. Such a model should contain 
the basic physics involved: Each magnetic nucleus is sen- 
sitive to local magnetic polarization of electrons which 
changes the effective field. This manifcstsfl as different 
(Knight) shifts in the NMR absorption frequencies lead- 
ing to inhomogeneous line broadening and different relax- 
ation times. While it is obvious that a close connection 
exists between the output of an NMR experiment and the 
geometry and disorder of the particle on which the elec- 
tronic properties depend, we are still far from a universal 



quantitative description of the effect. 

In this paper we address the problem of how the size of 
a particle contributes to the NMR lineshape. Because of 
quantum interferences, the de^uls of shape and disorder 
of the nanoparticles manifesioi3 as a Wigner-Dyson or 
Poisson distribution of electronic energy levels depend- 
ing on whether the states are extended (metallic phase) 
or localized (dielectric phase). In both cases, there are 
strong fluctuations in the local density of states (LDOS) 
and hence in the local Pauli susceptibility. Since samples 
used in earlier experiments had a distribution of parti- 
cle sizes, the inhomogeneous line shapes were assigned 
to this distribution. Different sizes would not only give 
a different number of electrons at the Fermi energy (av- 
erage density of states) but also involve a fluctuation in 
the parity of total electron number giving a further con- 
tribution to the inhomogeneity. However, these efforts, 
based in a canonical ensemble, do not give a consistent 
explanation for the observed lineshape valid for all the 
systems studied. A recent break-through was proposedllll 
with the use of grand canonical ensemble which allows 
for fluctuations in the particle charge. This picture, 
which is valid in many experimental situations, empha- 
sizes the fact that inhomogeneities manifest already in 
a single particle. Therefore inhomogeneous broaden- 
ing occurs jfiven when particles are exactly equal. This 
is observedCI for the ligand-stabilized metal-cluster com- 
pounds PtaogPheuggOao. Estimations of this shape were 
performed using an electronic structure represented by a 
non-linear sigma modeltll and by a Random Matrix The- 
ory (RMT)E£l. Both models produce the same results for 
the LDOS fluctuations. One could visualize the model 
Hamiltonian implicit in the RMT theory, considering a 
set on M "atomic orbitals" whose interactions, not con- 
strained to a given lattice, are described by the Hamilto- 
nian matrix elements Ti-ij, which are random and satisfy 
a gaussian distribution ' P(H) = C exp[-TrHV(2cr^)], 
with C a normalization constant and ^ M to keep 
a finite band width B. Since the analytical solutions 
are obtained in the limit M — > oo, each "orbitaF sees an 
equivalent " mean field" environment. While solutions de- 
scribe fluctuations in the LDOS, they can not distinguish 
surface from bulk "orbitals" or the explicit geometry of 
the particles neither describe the transition to the local- 
ized regime. Hence, while these analytical models contain 
the essential physics of level statistics, they are not ex- 
pected to provide a detailed characterization of different 
materials and particle geometries. In order to overcome 
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this limitation, in this work we use a very simple tight 
binding model for the electronic structure that is a first 
step toward a more precise description of the nanoparti- 
cles in terms of Linear Combination of Atomic-Orbitals 
(LCAO). A similar model was used recentlyB13 to de- 
scribe the cross over between Wigner-Dyson and Poisson 
statistics. We will also show that resorting to appropri- 
ate calculation tools (e.g. Matrix Continued Fractions) 
such models can be used in practical calculations. 

II. THE MODEL 

The electronic system is described by the Hamiltonian 

M MM 

n=Y, E.c+c, + ^ ^(lS-,^c+c, + y.,,c+c,.) (1) 

i—l j>i 2—1 

where Ei is the energy of the s-like atomic state cen- 
tered at nuclei position orthogonalized atomic orbital 
ipi{r) = ip{r — Ti), and Vj^i are real numbers accounting 
for the kinetic energy of the electrons. There are AI — 
il/a^ orbitals in the volume O = L1L2L3. To model the 
disorder we assume that Ei^s are uniformly distributed 
around the mean value Es in the range [—W/2,W/2] 
(Anderson model). In a Hartree mean field description, 
charge transfer energies are accounted for conserving the 
structure of (P. The presence of the magnetic field H 
leads to the unitary symmetry and can be included in 
the Hamiltonian through the Peierls substitution. How- 
ever, for experimentally accessible fields {H < lO^G) and 
typical particle sizes (L > lO^^cm) its effect is negligi- 
ble {HL'^e/ (ch) < 1 with e/ (ch) = 2.4 x lO^G-icm-^) 
because degeneracies are already broken by the disorder 
and we restrict our calculation to orthogonal symmetry. 
In principle, interactions with other particles produce a 
decay of the eigenstates and enter through an inhomoge- 
neous bjxiadening of the ionization energy of orbitals at 
surfacealj {Ei — > Ei — iF^), accounting for electron tun- 
neling toward neighbor particles. While this is not much 
complication in a numerical treatment of the problem, 
to allow for comparison with previous works, we adopt 
here a simpler description in which F^ = rjo at every site. 
Properties of the single particle excitation spectrum are 
contained in the retarded (advanced) Green's function 

^^^-\^e+[{-~)i^-E,y 

where V'fc(r) = o.k(i)fi('>^) and Ek are the exact eigen- 
functions (molecular orbitals) and eigenenergies for the 
isolated particle. 77 is a natural broadening of the elec- 
tronic states and the (— ) sign corresponds to the retarded 
Green's function. The local densities of states per atom 
(LDOS) at «-th site is evaluated as: 

iV.(£) = -(27ri)-MG5(£)-G^,(e)], (3) 



from which the relevant contribution to the density of 
states per unit volume at the «-th nucleus, iV(£, ri),is 
obtained. In presence of an external magnetic field, elec- 
trons polarize proportionally to a local Pauli susceptibil- 
ity 

Ar(e,r,)[--]de, (4) 

with f(e) is the Fermi function. Numerical integration of 
(^) for finite temperatures can be very time consuming 
due to the singular nature of the spectrum which would 
require a LDOS calculation for each energy. In order to 
avoid integration, the finite temperature effects can be 
taken into account by taking 

_ r9f ]_ arkBT 

^de^^ TT{e-eFY + {aTkBTf' ^' 

which in practice becomes equivalent to put ksT ~ 
in Eq.^ and use an energy uncertainty for the electronic 
levels 

r] ^ r]o + ksT, (6) 

where the coefficient = 1 gives a good approximation 
approximation for Eq. (||). The Knight shift at a given 
nucleus is 

Ac^(r,) = J—^p{r,)H, (7) 

where J is the hyperfine constant and g/i the electron 
magnetic moment. Then, Auj{ri) will follow a probability 
distribution proportional to that of the electron LDOS. 
The NMR absorption line is inhomogeneously broadened 
by this distribution and is obtained by summing up all 
magnetic nuclei. The resonance line finally becomes 

liij) = AY^S[lu-luo~ Auj{r,)], (8) 

i 

where uq is the non-shifted position of the resonance in 
a dielectric material (as in a salt). That is, the observed 
NMR line shape is the sum of narrow lines with and 
inhomogeneously distribution in frequencies that follows 
that of the LDOS. In order to allow comparisons between 
systems, it is convenient to define the adimensional fre- 
quency variable x — (lu — uJo)/ {'^k — ^o), where ojk is the 
peak resonance in the bulk metal. This gives x = 1 for a 
metal and a; = for the dielectric phase. 

The intensity can be normalized: ^ I{x)dx = 1. A 
quantity of interest to characterize the line is the Knight 
shift Xm defined as I{xm) = max[/(x)]; another is the 
first moment xi = J xl{x) dec, whose difference from Xm 
is a measure of the asymmetry of the line. Notice that if 
we are working within a grand canonical ensemble in an 
infinite system, we expect that even when an electronic 
density at some nuclei is diminished by interferences this 
is compensated by an increase in others. In the analytic 
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modelS this property results in = 1. However in a 
MO model of the cluster this is not necessarily true. In 
a cluster, the ensemble average DOS is not a monotonic 
function of energy but an oscillating function with local 
maxima at the eigenenergies of the ordered cluster. This 
dependence on on the cluster size is an easy analytical 
resultllil for a Lorcntzian distribution of site energies. 



III. MATRIX CONTINUED FRACTION 

The exact solution of Hamiltonian (|l]) is a formidable 
task. Rough features of|-t|he LDOS can be obtained 
by the Recursion Methodic. This involves the obten- 
tion of a tridiagonal basis for (|l|) and the evaluation of 
the local Green's function (||) as a continued fraction. 
When used within the RMT this method allows analyti- 
cal solutionaij. We resort to it for the evaluation of the 
density of states per orbital in an infinite ordered system, 
NQ{eF,Ti)- From this quantity the typical level separation 
for a cluster with M orbitals is calculated: 



A = l/[M ■No{eF,v)]- 



(9) 



However, for small rj (low temperature), it is difficult 
to obtain the precise fluctuations in the LDOS for fi. nitO | 
clusters, because the intrinsic numerical instabilitiealli 



of the tridiagonalization procedure distort the fine de- 
tails of the LDOSEj. It is therefore compelling to resort 
to a more robust method to describe the infinitesimal 
details of the density of states. This is the Matrix Con- 
tinued Fractions (MCF) methodEHl for the calculation of 
the local Green's function. Its basic idea is to exploit the 
short range interactions in Hamiltonian ([^) by indexing 
states in a way that subspaces representing layers inter- 
act through nearest neighbor subspaces. In matrix form: 
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, (10) 



where O's are null matrices, E's in the diagonal repre- 
sents intra-layer interactions while the only non-zero off- 
diagonal matrices and V„^„-|-i connect nearest neighbor 
layers. Detailed structure of the submatrices depends on 
the lattice, for the cubic structure V„_„±i = VI, with 
1 the identity matrix. The local retarded Green's func- 
tions connecting sites i and j within the n-th layers are 
arranged in a matrix 



G^,n (e)= [(e + i??)l-E„ 



(11) 



where the matrix self energies and are cal- 

culated in terms of Matrix Continued Fractions (MCF) 
defined through the recurrence relations: 



n,n±l 



(e + ir;)l-E„±i,„±i -S^±± 



(12) 



which are calculated with the boundary conditions: 
= ST = 0. The intrinsic stability of the method 
manifestsEJ in the fact that a modification of a layer at 
distance L away from a given one will produce an ex- 
ponentially small variation of the local density of states, 
6N ~ exp[— 7iL] . Here 71 is the minimum Lyapunov 
Characteristic Exponent. In the localized regime {W > 
Wc, rj = 0) 1/71 coincides with the localization length 
while in the ordered case. [W <C Wc, 77 > 0) it gives 
a phase breaking lengthEj L^{ri) which washes out the 
finite size effects when L > La,. 



IV. RESULTS 

In this work we consider a simple cubic structure form- 
ing clusters of cubic shape with sizes L ranging from 5 
to 15 layers. That makes clusters of up to 3375 orbitals. 
We chose an arbitrary Fermi energy ep — Es + 4:\V\ 
which is far enough from the van Hove singularities for 
that lattice and hence allows a comparison with results 
from the non-linear sigma model which implies nearly 
free electrons. However, with —V w leV = 1.16 x 10'*K, 
it gives a reasonable scale to describe an s metal such as 
Cu. Disorder is small enough {W = 0.5 \ V\) to be on the 
metallic side of the metal insulator transition but strong 
enough to provide randomness of the energy levels. Fur- 
ther sampling is achieved by considering an ensemble of 
10 Fermi energies in a range Sep = 0.5 \ V\. This mini- 
mizes the effect of finite size correlations in the electron 
energy spectrum. 

Figure 1 shows the normalized distribution of LDOS 
for a cluster of size L=15 with three temperatures pro- 
ducing effective broadenings rj/V = 0.05, 0.0067 and 
0.001; which according to Eq. (^ can be interpreted 
as a sequence of decreasing temperatures. The number 
of occurrences I{x) of each normalized density of states 



(13) 



is proportional to the NMR absorption. Each member 
of the ensemble, containing four disorder configurations, 
is normalized according to the corresponding Fermi en- 
ergy and line broadening, a dependence writen explicitly 
in Eq.(p^. The same results were obtained with a sin- 
gle cluster, since the system is big enough to provide a 
configurational average. It is clear that by lowering the 
temperature the Knight shift decreases from its metallic 
limit Xm = 1 toward Xm — corresponding to a salt. 
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Simultaneously the line becomes increasingly asymmet- 
ric j-jXhese results are consistent with the experiments in 
CuEB and Pti. 

In Fig. 2 we show the variation of the lineshape with 
cluster sizes (L=10, 8 and 6) for a fixed natural broad- 
ening of rj/V — 0.018. Again it confirms the general 
trend experimentally observed in Cu and the predictions 
of the analytical models, i.e. the band tends to broadens 
and maxima shift to lower frequencies as particle size or 
temperature is decreased. 

In order to capture the universal nature of the com- 
petition between cluster size and temperature suggested 
by Figs. (1) and (2), we plot in Fig. (3) the value of 
the adimensional Knight Shift as a function of the scal- 
ing parameter 77/A = a/(27r) which measures the ratio 
between the level broadening and the typical level 
aration. a is the scaling parameter defined in ref. 
For completeness we also include in the plot a continuous 
curve obtained from the analytical lineshape obtained in 
that work for the unitary ensemble. As expected for a 
given 77/ A their results tend to give a higher Xm than 
our numerical values. This is consistent with the fact 
that the breaking of time reversal symmetry produces a 
stronger level repulsion diminishing the fluctuations. 

However, there is no universal scaling function that fits 
the complete lineshape for all systems. This is because in 
our model there are surfaces. Although we did not assign 
different ionization energies to surface orbitals (which can 
produce Tamm states), local densities of states at sur- 
face region (about 6/L of the total) can differ substan- 
tially from that at inner atom^J. Sites at the surfaces 
are connected to fewer neighbors and this manifest on 
a decrease in the wings of thc-LDOS as a consequence 
of lower gain in kinetic energycj, which is compensated 
with an increase of the LDOS at energies around Es- 
It should be emphasized, however, that the sign of the 
change in the LDOS at surfaces depends both, on the 
lattice structure and on the value oi ep', e.g. if = Es 
the surface LDOS would be higher than its value in the 
bulk. In any case, this produces an additional fluctu- 
ation in the LDOS which is more important for smaller 
systems. These effects attenuate within a few lattice con- 
stants and ape, consistent with the exponential healing 
length modeled. In Fig. (4) we separate the contribu- 
tions to the distribution of LDOS (linewidth) from sur- 
face and inner sites in a small cluster. Similar separations 
of maxima are obtained at higher 77/ A. To emphasize the 
importance of this effect on the line broadening in small 
particles, we show in Fig. (5) the scaling function satis- 
fied by function of the adimensional parameter 
ry/A for a 5 X 5 X 5 cluster in which we also discriminate 
the contributions of surface and inner sites. It is clear 
a general tendency toward lower values of the Knight 
shift for nuclei at the surfaces. Maxima do not coincide 
even when 77/ A is increased. In small clusters this effect 
manifests as an intrinsic linewidth persisting at high tem- 
peratures, a tendency that departs from the scaling lajw. 
This is the situation observed in recent experimentE 



on Ptaog and NisgPtg particles, where linewidths do not 
shrink with increasing temperature. Our results suggest 
that size and temperature dependence of line shapes can 
be used in jaonjunction with other NMR techniques such 
as SEDORH to allow the experimental identification of 
electronic densities at surface regions. 



V. FINAL REMARKS 



In agreement with analytical theoriesllillj for meso- 
scopic fluctuation, our results show that distribution of 
LDOS (and hence NMR line) shifts and asymmetrically 
broadens as temperature decreases and particle size di- 
minishes and can be described by a universal scaling law 
which is consistent general trends observed in Cu par- 
ticles and the size dependenceciZl of Ptaog and NiagPte 



particles, although the origin of the disorder in the last 
cases is not quite clear. In this case a finite broadening is 
obtained even for high temperatures. These departures 
from scaling theory in very small particles, is understood 
in our model as due to the presence of surfaces, a feature 
that the analytical theories do not contain. Therefore, 
our LCAO model contains the main features of meso- 
scopic fiuctuations and models the exponential healing 
of surface effects. Due to the simplicity and speed of the 
MCF calculation, it allows for further improvements of 
the model Hamiltonian which should lead to a more de- 
tailed evaluation of the electronic structure. Also orbital 
susceptibility could be quantified within these models. 
Since Cu has also a contribution to the linewidth due to 
quadrupolar interactions, and Pt, having a d-band con- 
tribution to the Fermi surface, has an important core po- 
larization susceptibility, they are far from ideal systems. 
A good candidate for a full quantification the described 
behavior is Ag, with spin 1/2 nucleus and s electrons. 
Therefore, further development along these lines would 
allow for: 1) a better quantification of NMR experiments; 
2) the investigation of the electron states determining 
the chemical kinetics at surfaces and 3) the progress in 
the understanding of basic phenomena at a mesoscopic 
scale. While the first issue was the main one addressed in 
this article we would like to add a few comments on the 
other two. Since in a frontier-orbital picture of bonding 
to metal surfaces, the states around the Fermi energy are 
those that determine the chemical reactions, it is of great 
interest to determine how these are modified by chang- 
ing the size and disorder of the catalyst particlen3. In 
particular a metal-insulator transition in the low dimen- 
sional (Tamm) surface band could play a very relevant 
role in chemical kinetics. Therefore, the theoretical and 
experimental methods described in this paper can be con- 
veniently adapted to suit this investigation. Finally, the 
interference phenomena lying in the core of mesoscopic 
physics could be controlled to obtain a fine tuning of 
many static and dynamical properties of materials. In 
particular, the present technique could be used to com- 
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plement current studieai3 that attempt to quantify the 
interplay of disorder, magnetic field and temperature in 
the metal-insulator transition. 
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VII. FIGURE CAPTIONS 



Figure 1. Typical distribution of occurrences 
/ of LDOS (or NMR lineshape) as function of the 
adimensional frequency x = {uj ~ ujq)/{ujk — ^^o) = 
N{eF,ri)/^o{£F,''l) evaluated for an ensemble of 4 clus- 
ters ofl5xl5xl5 orbitals sampling 10 energies around 
ep — Es + 4:\V\ at different temperatures: rj/V = 
ri^ + keT = 0.05, 0.0067 and 0.001 The disorder is 
W = 0.5 

Figure 2. Variation of NMR lineshape / for an en- 
semble of 10 particles of sizes L — 6, 8 and 10 lattice 
units, with 10 energies each. Temperature is such that 
1] — 0.018 \V\ , and other parameters are as in Fig. (1). 

Figure 3. Scaling curve for the Knight shift Xm for 
particles of different sizes as function of the adimensional 
parameter rj/A. The continuous curve is obtained from 
the analytical solution for a unitary ensemble of ref.EJ. 

Figure 4. Surface and bulk contribution to the line 
shape in a small cluster (M = 5x5x5), calculated with 
an ensemble of 20 different disorder configuration and 10 
energies each. Temperature is such that rj/V = 0.031. 
The dotted plot corresponds to normalized contributions 
of orbitals in the surface of the cluster. Full line plot is 
from inner sites. 

Figure 5. Surface and Bulk contribution to the 
Knight shift (or maximum in the line) is shown here for 
a small cluster with varying rj/V = rjo + ksT evaluated 
from curves as those in Fig. (4). Open triangles represent 
inner nuclei while open circles arc those on the surface. 
Filled squares represent the Knight shift evaluated from 
all sites. 
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